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We investigate the production of gamma-rays in the inner regions of SS 433 and in its 
interaction between its jets and the W50 nebula. We estimate the VHE emission that can 
be generated within the jets at distances > 10 13 cm from the compact object. We also 
apply a theoretical model of the jet/medium interaction to SS 433/W50. We compare the 
predicted fluxes to observations of SS 433 at TeV energies, and derive new constraints 
of the physical properties of this system. 
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1. Introduction 

Relativistic jets from a stel lar c ompact object were discovered for the first time 
in the microquasar SS 433 The system is composed of a 9 Mq, black hole 
orbiting a 30 M Q A3-7 supergiant star with orbital radius ~ 79 Rq and period 
■Porb ~ 13.1 d. Located at a distance of 5.5 ± 0.2 kpc^, the system shows relativistic 
jets with a velocity of 0.26c^. The jets precess with a period P prc ~ 162 d in cones 
of half opening angle 9 « 21° and are inclined by i ss 78° with respect to the 
observer line of sightM SS 433 is the only X-ray binary system in which hadrons 
have been found in the jet. Clouds of plasma with baryonic content are observed at 
large distances, suggesting that atomic reheating is working at zj ot < 10 17 cm from 
the compact object. Furthermore, a continuous regime of supercritical accretion 
onto the black hole is accomplished in SS 433. This could explain the large kinetic 
power thought to be transported by the jets, ~ 10 39 erg s _1 (see, e.g. Ref. [8}. 
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Fig. 1. Radio map of the extended W50 radio nebula (blue contours, extracted from Dubner et 
al. 1998) and the different regions where particle acceleration and emission could take place. Top: 
sketch of the SS 433 central engine. The jet acceleration region Az where the optical depth is < 1 
is represented by gray-shadowed cones. See text for details. Bottom: East and West interaction 
regions. The forward, reverse and reconfinement strong shocks are displayed. See Bordas et al. 
(2009) for a complete description of the formation of the interaction structures. 
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2. VHE observations of SS 433/W50 

SS 433/W50 has been extensively observed by the HEGIuP, MAGIcPl anc j 
CANGAROO-lfn] Cherenkov telescopes. No gamma-ray signal has been found, 
and upper limits at different energy thresholds have been established both for the 
inner system and the different interaction regions (the el, e2, and e3 regions at the 
east and the wl, w2, pi, p2, and p3 at the west; see Refs. 19) HU1 and ITT1 for details). 
In the table below we list some of the reported upper limits. MAGIC observations, 
in particular, took into account the strong absorption due to the periodic compan- 
ion eclipses as well as the attenuation due to the precession of the accretion disk 
envelope (see Ref. [12] for a detailed study) . 



Table 1. Observations of SS 433 and the jet/W50 interactions in the TeV regime by the MAGIC, H.E.S.S. and CANGAROO-II 
Cherenkov telescopes. The total observation time, the energy threshold E t ^ and the reported upper limits to the integral flux 
(f> u \ (in units of 10 — 12 ph cm -2 s — 1 ) are listed. 



Source 


IACT 


obs. time (h) 


Eth (GeV) 


4>ui 


Source 


IACT 


obs. time (h) 


E th (GeV) 


<t>ul 


SS 433 


HEGRA 


96.3 


800 


8.93 


Pi 


CANGAROO-II 


85.2 


850 


1.5 




MAGIC 


6.6 


250 


4.30 


p2 


CANGAROO-II 


85.2 


850 


1.3 


el 


HEGRA 


72.0 


800 


6.18 


P3 


CANGAROO-II 


85.2 


850 


0.79 


e2 


HEGRA 


73.1 


800 


9.18 


wl 


HEGRA 


104.9 


800 


6.65 


e3 


HEGRA 


68.8 


800 


8.96 


w2 


HEGRA 


100.7 


800 


9.00 



3. Gamma-ray production in SS 433/W50 
3.1. Inner regions 

We have estimated the Inverse Compton (IC) VHE fluxes produced at the borders 
of the binary system (see the upper inset in Fig. 1). At distances Zj et < 10 13 cm 
along the jet, the opacity to gamma-ray propagation due to the companion and 
accretion disk photon fields should make the optical depth r 7 > 1. We consider 
therefore acceleration/emission regions beginning further away than that point. The 
free parameters in our model are the non-thermal acceleration fraction q re \ and the 
accelerator/emitter size Az — 2j C t,f — z jct,i- Table 2 lists some of the parameters 
that remain fixed in our model and their respective values (parameters from the 
interaction model are also displayed, see below). 

We assume the dominance of the IC channel over any other relevant energy 
loss mechanisms. To estimate the 7-ray fluxes we approximate both the companion 
star and the disk envelope as point-like sources of iso tropic black-body radiation 
at temperatures T=8500 K and 21.000 K, respectively I 13 | 14 | yye use ^ ne iq cross 
section including both Thomson and Klein-Nishina regimes for this process, since 
7c ^ 7kn = (4eo) _1 , where 7 e is the electron Lorentz factor andeo xm cC 2 ~ 2.7 K^T 
is the peak energy of each corresponding photon field. We take a given fraction of 
the total jet power to be delivered to a leptonic accelerated plasma that follows a 
power-law distribution, N e oc 7,T P , with a spectral index p — 2. A maximum particle 
Lorentz factor ^™ ax = 10 6 is used in our computations. 
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3.2. Interaction regions 

Three strong shocks are formed when the jet impacts with the surrounding nebula 
(see Fig. 1, bottom). A forward shock, propagating into the medium and enclosing 
the shell region; a reverse shock directed into the jet material and inflating the 
cocoon, and a recollimation shock that appears when the pressure in the conical jet 
equals that of the cocoon. We have studied the non-thermal emission coming from 
these three regions through synchrotron, IC and relativistic bremsstrahlung (the 
last only in the shell, where it is relevant) emission channels. The model assumes a 
self-similar growing of the interaction structures to estimate the physical properties 
of each shocked region. The reader is referred to the work presented in Ref. [15] for a 
complete model description. Some relevant model parameters are listed in Table 2. 



Table 2. List of the parameters that remain with a constant 
value in the analytical model. 



Parameter 


Symbol 


Value 


Jet kinetic power (erg s — 1 ) 


Qjct 


10 39 


ISM density (cm -3 ) 


"ISM 


1 


Source age (yr) 


*MQ 


5 x 10 4 


Jet Lorentz factor 


Tjet 


1.04 


Jet semi-opening angle (°) 




1.2 


Luminosity companion star (erg s — 1 ) 




10 39 


Self-Similar parameter 


R 


3 


Magnetic cquipartition fraction 


V 


0.1 



4. Model results 

Different integral luminosities from the inner system regions for energies £L > 
250 GeV are obtained if gamma-rays are generated at distances Zj e t,i ~ 10 13 ~ 14 cm 
up to Zj e t,f ~ 10 17 cm. We note that taking Zj e t,f much higher than 10 17 cm does 
not significantly change our results, since at these distances the photon field en- 
ergy density is negligible (we do not consider the contribution from the CMB in 
our calculations). The integral fluxes produced along the accelerator/emitter region 
/\ z accei ag a f unc ti on f th e injection point and the acceleration efficiency are shown 
in Fig. 2. The MAGIC upper limit < 4.3 x 10~ 12 ph cm -2 s^ 1 for E 1 > 250 
GeV for the central regions of SS 433 is also displayed. 

Concerning the interaction of SS 433 with the W50 nebula, we show in Fig. 3 the 
SED computed for the shell, the cocoon and the reconfinement interaction regions. 
Relativistic bremsstrahlung is the most efficient gamma-ray production mechanism 
in the shell zone, reaching bolometric luminosities at a level ~ 10 33 erg s _1 . IC 
is the main radiation channel at high and VHE in the cocoon and reconfinement 
regions, with bolometric luminosities at the level of a few x 10 31 erg s _1 . The integral 
gamma-ray flux for energies E > 800 GeV computed from the SED is found to be 
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Fig. 2. MAGIC integral flux upper limit at E 7 > 250 GeV (horizontal blue line), together with 
the total flux obtained through IC process for three different values of the acceleration efficiency 
q IC \ and acceleration injection height Zj c t,i ( we recall that leptons are accelerated along Az?™ 6 ' = 
Zjct.f — Zjet.i, and we have used Zj e t f = 10 17 cm in all cases; see text for details). The region at 
the left side of the red-dotted line has optical depth values > 1 . 




Fig. 3. SED of the interaction regions obtained for the parameters listed in Table 2. The sum 
of the contributions from the shell region, the cocoon and the reconfinement region is displayed. 
Relativistic bremsstrahlung is important only in the shell region, due to the low particle densities 
in the jet and cocoon. The contribution of only one jet impacting on the nebula is showed. 



slightly lower than the upper limits for the different interaction regions listed in 
Table 2. 
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5. Discussion 

We explored the 7-ray emission produced by relativistic electrons both in the inner 
jet regions and in the jet/W50 interaction regions. The gamma-ray fluxes obtained 
for the central regions depend linearly on q rc i and the acceleration region within the 
jet, which are nonetheless poorly constrained under a theoretical point of view. A 
value of q le \ far larger than <~ 1CP 6 for Az^ t cel starting at Zj ot ,i ~ 10 13 cm is ruled 
out within our flux estimations and the reported upper limits, while q Te \ > 10~ 5 is 
not allowed provided that the acceleration region Az?^ t cel starts at Zj C t,i > 10 14 cm. 
Acceleration processes could be much less efficient than expected, though it seems 
unlikely since acceleration is actually required to explain the non-thermal emission 
at lower energies. In addition to the IC process, relativistic bremsstrahlung and 
SSC could be important at the very inner regions near the jet base, where magnetic 
fields and ion densities are the highest. However, the opacity in these zones should 
make their contribution at VHE difficult to observe. 

No 7-ray signal has been found from the SS 433/W50 interaction regions. Nev- 
ertheless, the huge jet kinetic power and the high density of the W50 nebula, make 
this source one of the best microquasar/ISM interaction TeV-emitting candidates to 
be pursued. Our analytical treatment predicts fluxes that are slightly below the cur- 
rent telescopes sensitivity. Nevertheless, the improvement of the Chcrenkov facilities 
(mostly the upcoming MAGIC and H.E.S.S. upgrades) could bring the opportunity 
to detect such interaction structures in the TeV regime for the first time. 
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